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\Q ; ABSTRACT 

^J , We explore the benefits of using a passively evolving population of galaxies to measure 

^— V ■ the evolution of the rate of structure growth between z — 0.25 and z = 0.65 by 

^vq I combining data from the SDSS-I/II and SDSS-III surveys. The large-scale linear bias 

of a population of dynamically passive galaxies, which we select from both surveys, 

is easily modelled. Knowing the bias evolution breaks degeneracies inherent to other 

methodologies, and decreases the uncertainty in measurements of the rate of structure 

k^ . growth and the normalization of the galaxy power-spectrum by up to a factor of two. If 

^ ■ we translate our measurements into a constraint on as{z — 0) assuming a concordance 

C^ I cosmological model and General Relativity (GR), we find that using a bias model 

improves our uncertainty by a factor of nearly 1.5. Our results are consistent with a 
flat A Cold Dark Matter model and with GR. 
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1 INTRODUCTION | 2010l : iReid et al.ll2010l : iBlake et al.ll2011al lbll3: IConlev et al.l 

I2OIII ). The physical reason behind such an acceleration re- 
Current observational evidence points towards a Universe . , 1 j_ 1 1 ^ j_- 1 1 i- 

mams poorly understood, and potential explanations range 
that IS undergoing an accelerated expansion (see e.g. „ . , 1 • 1 ^ j_ 1 -^ ^ 

i Kessleretal l & Umanullah et al I lioiol- I Percival et al l "" '™P cosmological constant or vacuum density, to 
i-Lvctaoicr cL ai.i i-^uuci . i-Tiiiiaiiuiiaii cl ai.i izuin . ir^cicivai cl ai.i i.ni • . t ^ -i tt- 
modmed gravity models or an mhomogeneous Universe cre- 
ating the illusion of an acceleration. Distinguishing between 
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such physical explanations is a key goal of modern Cosmol- 
ogy- 

Redshift-space distortions (RSDs) are a key observa- 
tional tool for understanding Dark Energy as they trace the 
matter velocity field via the peculiar velocities of galaxies. 
They allow a measurement of the growth rate of structure 
via an enhancemen t of the clustering power along the line of 
sight (jKaiserll 19871 1 . RSDs are powerful discriminants of dif- 
ferent physical models for Dark Energy, as models that share 
the same expansion hist ory often predict differe nt growth 
rates of structure, / (e.g. iLinder fc Jenkinsll2003l ). 

Large-scale clustering measurements yield a direct mea- 
surement of fas and has, where / is the logarithmic deriva- 
tive of the linear growth factor D{z) with the scale factor, 
/ = d log _D(2;)/d log a, as is the variance of the matter den- 
sity field at a scale of 8 h~^ Mpc, and b is the large-scale 
linear galaxy bias. These results must be coupled with in- 
dependent measurements of b or as to yield an estimate 
of the growth rate, which often requires further assump- 
tions: galaxy bias measurements are notoriously difficult, 
and measurements of erg often need to be extrapolated in 
redshift. Higher-order clustering measurements can also be 
used to break the degeneracy between g alaxy bias and cos- 
molo gy (see e.g. Bernardeau et al. 2002l : IZheng fc Weinberg 
|2007|) which has been investigated with galaxy data (see 
e.g. 'Pan fc Szapudi"2005'; 'Caztafiaga et al.ll2005l : IRoss et all 
[2008; Marin 2011; McBridc et al. 20ll|). Obtaining precise 
constraints from higher order moments is challenging, and 
this work serves as a complement to such investigations. 

In this paper we explore the gain if one knows the ex- 
pected evolution of the bias for a sample of galaxies. For 
a passively evolving sample (i.e., no merging) t he bias evo - 
lution ca n be computed using the formalism of iFrvl (| 19961 ) 
(see also lTegmark fc Peebleslll998l ). This formalism models 
a population of galaxies as being formed by a non-linear 
process at some time in the past, and subsequently evolving 
with the velocity flows set up by the matter density field. 
To first order (i.e. in the linear regime), a simple model for 
the evolution of bias can be constructed. This formalism as- 
sumes the continuity equation for the galaxy density field, 
which conserves the number of galaxies as a function of red- 
shift, and imposes the need to select a dynamically passive 
sample. 

We obtain a passivel y evolving sample o f galaxies via 
the method described in iToieiro et al.l (|2012l ). which pro- 



vides weights and carefully matched galaxy samples span- 
ning SDSS-I/II and III. Galaxies are weighted according to 
the volume in which they would be visible across the two 
surveys, and this matches the samples from an evolutionary 
point of view - SDSS-III galaxies seen through most of the 
SDSS-I/II volume are more likely to be their progenitors, 
and SDSS-I/II galaxies seen through most of the SDSS-III 
volume are more likely to be the evolutionary products of 
SDSS-III galaxies. One can then assess the consistency of 
this weighted sample with a dynamically passive model by 
computing the evolution of the number and luminosity den- 
sities - in a purely passive model, these should be constant 
with redshift. The most robust es timate of the merge r rate 
in the weighted galaxy sample of iToieiro et al.l (|2012f ) over 
0.2 < z < 0.7 yields a modest value of 2% ± 1.5% Gyr"\ 
establishing its suitability for our present study. 

When computing the large-scale clustering amplitude 



we weight each galaxy by its luminosity, and we construct 
samples at each redshift to have the same weighted lumi- 
nosity density. The luminosity weighting gives a large-scale 
power estimator that is less sensitive to galaxies within 
the sample merging between any two redshifts: i.e. merging 
events between galaxies in the same halo do not affect the 
relative contribution of the halos within which they reside 
to the overall bias of the sample, provided total luminos- 
ity is conserved in such a merger. This is only strictly true 
in the case of no loss of light to the intra-cluster medium. 
Nonetheless, weighting by luminosity will almost always be 
better than any weighting scheme that depends on the num- 
ber of objects in a halo - when two objects merge the relative 
contribution of a given halo to the overall clustering signal 
will be reduced by 1/2 if weighting by number. It follows 
that, provided that the overall loss of light is less than 50% 
of the combined light of the merging system, we have an esti- 
mation of the bias evolution that is less sensitive to merging 
of galaxies within the sample, and to which the Fry model is 
more applicable. The luminosity matching simply prevents 
selecting less luminous (and less biased) galaxies at different 
redshifts in case of merging. This would happen if one was 
to match samples on number density, for example. 

It is these careful matching and weight ing s cheme s 
that justify the us e of the bias evolution of iFrvl l|l996l ). 
IToieiro et al.l (|2012l ) further demonstrated that, assuming a 
A-Cold Dark Matte r (ACDM ) model and GR, the bias evo- 
lution of lFrvl ([1996) provides a formally good fit to the data. 
Whereas in itself such a consistency is no proof of either the 
cosmological model or of the bias evolution model, it is a 
result that confirms our interpretation of the evolution of 
the galaxies within the broad context of a firmly motivated 
cosmological model. In this paper we assume the expansion 
history and matter power spectrum of a flat ACDM universe, 
but we independently measure the growth rate of structure 
that gives the best fit to the data - which may be decou- 
pled from the energy density and need not follow GR. The 
added constraint from the bias evolution allows us to break 
the degeneracy between galaxy bias, growth rate and erg. 

Finally we benefit from working on large scales (30-200 
h^^ Mpc); the modelling of the matter power spectrum and 
RSDs on non-linear and quasi-linear scales is poorly under- 
stood and a further source of uncertainty fe.g. lReid fc White! 
I2OIH ). In this first analysis we ignore most non-linear effects, 
accepting that future extensions of this work (with larger 
samples of galaxies and better statistical errors) will require 
a more sophisticated treatment of such effects. Where re- 
quired we assume a flat ACDM cosmology with ilm = 0.25, 
and Hq = 70kms"^ Mpc"^ 



2 DATA 

The Baryon Oscillation Spectroscopic Survey (BOSS), 
as part of the Sloa n Digital Sky Survey (SDSS) III 
(jEisenstein et al.ll201ll ). increased the total SDSS-I/II imag- 
ing footprint to nearly 14, 500 sq. degrees; all of the imag- 
ing was re-process ed as part of SDSS Data Release 8 
lAiharaet al.l l201l'). In SDSS-I/II, Luminous Red Galax- 
ies (LRGs) were selected for spectroscopic follow-up accord- 
ing to the target algorithm described in lEisenstein et al.l 
(J200l|), designed to follow a passive stellar population in 
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colour and magnitude space. In SDSS-III, the BOSS target 
selection extends the SDSS-I/II algorithm to target fainter 
and bluer galaxies in order to achieve a galaxy number den- 
sity of 3 X 10"'' h^ Mpc~'^ and increase the redshift range 
out to z ~ 0.7. The spectroscopic footprint of the BOSS 
data used here covers 3275 sq. degrees of sky, and corre- 
sponds to the upcoming Data Release 9, which will mark 
the first spectroscopic data release of BOSS. A set of com- 
pr ehensive clustering an al yses of this samp le can be found 
in I Anderson et al.! ('2012"'): ^ Reid et al.l (1201 2'): Sanch ez et all 
(|2012h : lManeraet al.. (,2012 ') and lRoss et al.. (,2012 '). The tar- 
get selection algor ithms for the LRGs and BOSS are de- 
scribed in detail in lToieiro et al.l (|2012l ). BOSS target selec- 
tion consists of two separate algorithms - in this paper we use 
only the CMASS (for Constant MASS) sample, selected to 
be approximately stellar-mass limited, and targeting galax- 
ies mainly with z > 0.43 (Padmanabhan et al. in prep). 

"We split the data across four redshift slices: two slices 
of LRG galaxies centred eX z = 0.3, 0.4 and two slices of 
the CMASS galaxies centred aX z = 0.5,0.6 (Az = 0.1), 
with 44136, 30393, 39780 and 37883 objects respectively. At 
each redshift we select the brightest galaxies until a fixed 
luminosity density is reached. This corresponds to roughly 
95% and 40% of the LRGs and CMASS samples respectively. 



3 THE MODEL 

We describe the rcdshift-space galaxy correlation function 
^(/i,r) as in [Hamilton. (^199^): 



can be com puted from g 
tegrals (s eelHamiltorJ 



e(M, 



eo(r)Po(M) + 6(r)P2(M) + i^{v)P^{pi) (1) 



where r is the comoving separation in hT^ Mpc and \i is the 
cosine of the angle between a galaxy pair and the line of 
sight. Pi are the Legendre polynomials with Po = 1, -P2 = 
(3/i^-l)/2andP4 = (35ai'*-30ai^-|-3)/8. ^n is the monopole 
of the correlation function; the excess of finding a pair of 
galaxies at given distance r averaged over pairs observed at 
all angles with respect to the line of sight. The quadropole, or 
^ = 2, contains the next order of information, by effectively 
comparing the power along and across lines of sight. Current 
measurements of the octopole, or ^ = 4, are too noisy to 
yield useful constraints and are not included in our model. 
We model the redshift evolution a nd the amplitud e of the 
monopole and of the quadruple as (|Hamiltonlll9921 ): 



io^T^z) 



h\z)^\i{z)h{z) + \f{z) 



ol{z)i^(v) (2) 



Ur,z) = ~ |/(,)fe(,) + l/(^) 



ol{z)i^{T) (3) 



with a%{z) = as{0)D{z)/D{0) where we set D(0) = 1, and 



, ^ using a set of well-defined in- 
I992I). In t his paper we use the ^0^2 ('^) 
models of ISamushia et al.l (|2012r ). with Q,m = 0.25. 

We describe the three-equation system above with 4 pa- 
rameters consisting of b{zo) and three nodes for (Js{z), which 
we model using a a quadratic polynomial. The nodes are at 
Znode ~ 0, 0.3 and 0.6; we find that changing these nodes 
within this range does not affect our results significantly. 



4 THE MEASUREMENTS 

We estimate the correlation function from the data, ^i{r), by 
means of the lLandv fc Szalavl (|l993l ) estimator. We use 130 
bins in r, logarithmically spaced between 1 and 200 h~^ Mpc, 
and 200 linear bins in /i, between and 1. We use a random 
catalogue with the same angular mask as the data catalogue, 
and with an n{z) matched to that of the data but with 10 
times the number density. The non-trivial survey geome- 
try imprints a non-uniform distribution of pairs in /i on the 
data. We correct for this effect as in lSamushia et al.l (|2012h . 
by weighting each galaxy pair such that the weighted dis- 
tribution of pairs in fj, corresponds to that expected in the 
absence of a survey ma sk. We correct for ang ular and red- 
shift completeness as in [Anderson et al.l (|2012r ). 

We weight each gal axy by its luminosi ty and Vmatch 
weight as described in iToieiro et al.l (|2012l '). The Knatch 
weight preferentially selects galaxies seen across both sur- 
veys and more likely to belong to the coeval population of 
galaxies we wish to consider, and the luminosity weighting 
results in an estimate of the large-scale power that is less 
sensitive to merging within the sample - see Section [T] To- 
gether these weights ensure the bias model of Equation (|4} 
is applicable to our sample. 

For each of the redshift slices we compute ^0,2(7'), and 
use a simple 2-dimensional x^ minimisation to find the 
best fitting scale-invariant amplitudes, Ao,2(z), by writing 
(,o,2{r,z) = Ao,2(z)go^2('^)- To ensure a stable inversion of 
the covariance matrix, and to increase our signal-to-noise 
in each bin, we re-bin ^0,2(2) to 11 bins between 30 and 
200 h~^ Mpc. Re-doing the analyses using scales between 50 
and 200 /i~^ Mpc significantly increases our overall errors, 
but does not change our conclusions. 

We estimate the errors and their covariance by using 
mock simulations. We use the LasDamas mocks (McBride et 
al. in prep) to construct 80 independent realisations of ^0,2 
for the first two redshift slices (we sub-sample each mock in 
order to reproduce the n{z) in each slice). F or the last two 
redshi ft slices, we use 600 PTHalo mocks of iManera et al.l 
(|2012f) . and follow the same procedure. We include the co- 
variance between the multipoles in our fits. The CMASS 
mocks assume a slightly different ACDM cosmology and are 
heavily subsampled to match the data n{z); we scale their 
mean correlation function to match the data and apply the 
same factor squared to the full covariance matrix. 



b{z) = [b{zo)~l]^ + l, 



(4) 



where equation ^ follows the modelling of iFrvl (| 19961 ) for 
evolution of the large-scale linear bias. ^0^2 hold the infor- 
mation on the shape of the matter correlation function, and 



5 RESULTS 

We adopt a Markov Chain Monte Carlo (MCMC) technique 
to sample the posterior distribution of the parameters in our 
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model, given the data. We set uniform priors on our param- 
eters as follows: 1 < 6(zo) < 3.5 and < as{znode) < 1-5. 
The marginalised likelihood distributions of all our param- 
eters have fallen to zero near these boundaries. We use a 
stationary proposal density function, with a shape similar 
to the marginalised likelihood distributions of each of our 
parameters, which we investigate with a set of preliminary 
chains. In each step of the chain we update one parameter at 
a time, randomly chosen and all with equality probability. 
Our final chains have an acceptance rate of r; 15%, and our 
results and la intervals are robust to changes in the choice 
of the proposal and starting point; different choices for the 
proposal simply lead to lower acceptance rates. We adopt 
the mean value of each marginalised distribution as being 
the best-fit value for a given parameter, and we take la er- 
rors from the standard deviation of the same distributions. 



5.1 Passive model 

Fig. [T] shows the marginalised likelihood distributions for 
the free parameters in our model: b^g and as{znodes) (first 
two panels), as well as for the derived parameters: f{z)as{z), 
b{z)as(z) and f{z). We choose to present the distributions of 
the derived parameters at the centre of the redshift slices we 
use to measure the correlation function, but note that these 
are not independent. The correlation factor between adja- 
cent measurements of f{z) is high, between 0.84 and 0.92, 
but between the two furthest measurements, a,t z — 0.3 and 
z = 0.6, it is lower (0.147). The correlations of f{z)as{z) are 
similar. We show the best-fit values and la confidence inter- 
vals in Table [T] under the header of passive model. The co- 
variance matrix for our fitted parameters is given in Table [5] 
- this is the parameter set and covariance matrix that should 
be used for estimating likelihood surfaces. Fig. [2] shows in 
red our measurements of f{z)as{z) as a function of redshift, 
compared to measurements from the literature. 



5.2 Free growth model 

To place the results from the previous Section into context, 
we fit fas and has independently in each of the redshift 
slices. We continue to use equations ^ and ((3]), but now 
drop the constraint on the bias evolution given by Q. We 
use an MCMC similar to the one described in Section [SI 
adapted to refiect the different parameters in this model, of 
which there are eight. The evolution of /erg can be seen in 
the blue points of Fig. [21 and we show the full set of results 
in Table [1] under the header of free growth. We see a loss 
in precision of up to a factor of two in the estimation of 
f{z)as{z) and 6(z)crg(2:), when compared to the constraints 
obtained using the passive model. Note that the measure- 
ments quoted under free growth in Table [l] at each redshift 
are now independent. 



5.3 Constraining power 

As it is difficult to judge the constraining power of correlated 
measurements, we undertake the following exercise. Assum- 
ing GR and ACDM , we assess how well as{z = 0) can be 
constrained, using each set of points in Fig. [2] When us- 
ing literature data, we assume the likelihood surfaces to be 



gaussian, and in the case of multiple measurements we as- 
sume them independent. In the case of the measurements de- 
rived in this Letter, we use the best-fit as{znodea) values and 
their covariance. We show the resulting constraints in Fig. [3] 
The constraints from the passive model are approximately 
1.5 times better than a free growth model, an d competi- 
tive w hen compared to state-of-the-art results of Reid et al.l 
(2012!) on the full CMASS sample, and JBlake et al.1 l|2011di ) 
with Wiggle Z. 



6 SUMMARY AND CONCLUSIONS 

We demonstrate for the first time how using a passive sam- 
ple of galaxies can enhance the accuracy of the measurement 
of the growth rate, via the added knowledge of the evolution 
of the large-scale galaxy bias. Our results are fully consistent 
with a fiat ACDM model and GR. When compared to fit- 
ting bag and fag independently at each redshift, we find an 
increase in precision of up to a factor of two. If we translate 
our ACDM measurements into a constraint on (T8(0), as- 
suming ACDM and GR, we find that a passive model gives 
ag,{0) — 0.79 ± 0.045 which is a nearly 1.5 times improve- 
ment on the results obtained using a free growth model, 
as{0) = 0.785 ± 0.065. Furthermore, these constraints are 
comparable with those obtained using the measurement of 
iReid et all (|2012l '). cr8(0) = O.rSbto.lfo, whilst only using 
- 40% of the BOSS CMASS galaxies (but adding SDSS- 
I/II). This technique offers great potential, and it will deliver 
highly competitive results as BOSS gathers more data. 

A smaller statistical error in the measurements will re- 
quire a more sophisticated modelling of non-linearities in the 
treatment of RSDs, as well as a potential extension of the 
bias evolution model to accommodate a sample of galaxies 
that will be increasingly less dynamically passive as we ex- 
tend this work in luminosity and/or redshift. The obvious 
caveat is that we need to provide a convincing case that a 
sample is well match ed to passive evoluti on. For our sample 
this was provided bv JToieiro et al.l (|2012l '). 

With the right dataset and modelling, it is straightfor- 
ward to extend this technique to higher redshift, and map 
the growth of structure over a larger fraction of the age of 
the Universe. 
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Figure 1. Black curves in all panels show the marginalised likelihood distributions of our fitted and derived parameters. The fitted 
parameters are fe^p (first panel) and c8(-2nodea) (second panel, with z^ade. = 0, 0.3 and 0.6 from right to left). The derived parameters are 
/(z), f{z)(Js{z) and b(z)crs{z). Vertical solid red lines show the best-fit values, and the vertical dot-dashed red lines the Icr confidence 
intervals. Top right two panels show the measured value of Ao^2{z) (black circles) and la errors - the red line shows the best fit model. 
Dashed blue lines throughout show predictions from ACDM and GR, using the best-fit values for the fitted parameters. GR is perfectly 
compatible with our measurements of the growth rate. 
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Table 1. Summary of the results in this letter. The passive model corresponds to the model described in Section |3] using the bias 
evolution for passive galaxies. The free-growth model corresponds to the model described in Section 15.21 
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Figure 3. Constraints on 0-8(2: = 0) from tlie data points in 
Fig. [21 assuming ACDM and GR. The vertical shaded bar shows 
the constraints place d by the joint data analysis in WMAP7 
IjKomatsu et al.ll201lh . The constraints from the passive model 
are approximately 1.5 tim es better than a free growth model, 
and competitive relative to lReid et al.l l|2012l ') on the full CMASS 
sample. On the left we show the dataset used for each measure- 
ment. 



Figure 2. Evolution of /erg as a function of redshift for the 
paasive model and f ree growth. The blac k data points are from : 
iBlake et all ll20 11dVPcr cival et al] (12004'), 'Tegmark et al.' {2 00di) 
and lGuzzo et al . (2008); as collected by Song & Pcrcival (20091). 
We also show measur ements from Samushia ct alj 1 12012^ and 
from lReid et al.l 1120121). For com plet eness wc also show the mea- 
surements of iDavis et al.l l|201ll 'l and iTurnbull et aU ll2012l') from 
peculi ar velocities at 2: = 0.02, as compiled bv lHudson fc Turnbulll 
||2012| V The smooth solid line shows the prediction of ACDM and 
GR, using a WMAP7 cosmology with 0-8(2 = 0) = 0.81. 
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Table 2. Covariance matrix for the fitted parameters recovered 
from the MCMC chain described in Section |5] 
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